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In this paper we study the formation of giant electrons-positron pair halos around the 
powerful high energy extragalactic sources. We investigate the dependence of radiation 
of pair halos, in particular the spectral and angular distributions on the energy spectrum 
of the primary gamma-rays, the redshift of the source, and the flux of the extragalactic 
background light. 
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1. Introduction 

Interactions of very hig h en ergy (VHE) gamma-rays with the diffuse Extragalactic 
Background light (EBL)P^ as well as the the possible effects caused by propagation 
of secondary (pair-produced) electrons in the intergalactic medium (see e.g. Ref. [3]), 
make the high energy gamma-rays unique carries of cosmological information about 
the epochs of formation of galaxies and their evolution in the past (see e.g. Ref. IUG3 
[7l [HI [9]) . The method is based on the distinct features in the spectra of high energy 
gamma-rays arriving from distant extragalactic objects caused by interactions with 
EBL. 

A serious obstacle in practical realization of this method is our poor knowledge 
concerning the primary (intrinsic) gamma-ray spectra produced in the source It 
is believed that the well coordinated observations of gamma-ray blazars in different 
energy bands should allow gamma-ray astronomers to identify the radiation mech- 
anisms, derive the principal model parameters, and thus reconstruct the intrinsic 
gamma-ray spectra based on spectral and temporal properties of blazars (see e.g. 
Ref. [TT|) . Even so, for reliable estimates of the intergalactic absorption effect, and 
ultimately for derivation of the flux and spectrum of EBL, one should take into 
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account possible integral absorption of gamma-rays. The importance of the internal 
photon-photon absorption in blazars has been discussed by many authors, in par- 
ticular just after the discovery by COS B of the first gamma-ray blazar - 30273^21. 
The complex behavior of internal photon-photon absorption may result in a rather 
" irregular" deformation of the spec trum of primary gamma-rays (from steepening 
or even sharp cutoffs to hardening}^. This, as well as the effects related to the 
redshift of blazr J^, challenge the feasibility of reliable derivation of EBL based on 
absorption features in the gamma-ray spectra of extragalactic objects. 

Moreover, strictly speaking, the intergalactic absorption features contain infor- 
mation about the product of the diffuse extragalactic background radiation density 
w r and the Hubble constant Hq. Fortunately, it is possible to decouple w r and Hq 
by studying the spectral and angular characteristics of VHE 7 radiation from hypo- 
thetical electron-positron pair halos surrounding powerful nonthermal extragalactic 
objects. These giant structures are unavoidably formed around any extragalactic 
VHE source due to development of pair cascades initiated by interactions of pri- 
mary multi-TeV photons with EBlP^. 

2. Basic processes of formation of pair halos 

For formation of electron-positron pair halos around extragalactic objects, the 
photon-photon pair production (PP) and inverse Compton scattering (IC) are 
two principal processes which initiate and support the electromagnetic cascade 
development!^. The mean free pathlengths and the energy spectra of secondary 
particles produced at these interactions are described below. 

2.1. Pair production 

For the pair production at a collision of a gamma-ray photon of energy e with a 
target photon of energy lu the following threshold conditions is required: 

2(1 -/3' 2 )- 1 =ew(l-cos0), (1) 

where /?' is a velocity of the electron (positron) e ± in the CM frame and 6 is an 
angle between the momenta of two photons in lab frame (below we will use the 
energy of pairs in units of the electron rest mass, m c c 2 ). 

The mean free path length App of a gamma-ray of energy e traveling through 
the isotropically-distributed photon gas of differential energy density n(u)) is 

/•OO 

App(e) = / a P pn(uj)duj, (2) 

J Uth 

where a is the cross-section averaged over directions of interacting photons: 

^ /■!— 2/so 

a P p = - J (1 - [i)<jppdn . (3) 
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Fig. 1. Top panel: The cross section of photon-photon pair production app (in units of the 
Thomson cross section cry) calculated using the exact expression (solid line) and the approximation 
given by Eq. J5] (dashed line) . Bottom panel: the ratio of these two presentations of the cross-section. 



Here so = ecu, fi = cos 9 and app is the total cross section of the process (see 
e.g. Ref. OH) 



opp = 4 ct t(1-/?'V 
16 
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where <tt is the Thomson cross section. For integration of Eq. ([3]) we use the following 
approximation^^, 



_ 3<tt 



so + ~ In s - i + ^- ) \ ln( v ---n +- \ --n - 1 ) 
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(5) 



which provides an accuracy better than 3% as demonstrated in FigQ] 
The resulting energy distribution is described as^^ 



32w 2 e 3 



4e 2 4u;(e - 7)7 

x In 8u;e 



(€-7)7 
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Lue J (e — 7) 2 7 2 



(6) 



(e - 7)7 

where 7 is the energy of the electron (positron) (in m e c 2 units) and e = £|w«e 
is the photon energy in the lab-frame. The range of 7 is 



1-- <7< 1(1 



(7) 



Figure [2] shows the distributions of electrons (positrons) for different sq. 
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Fig. 2. Differential spectra of electrons and positrons from photon-photon interactions for 
different values of so: 10 (solid line), 3 (dashed line), 2 (dot-dashed line) and 1.5 (dotted 
line) . 



2.2. Inverse Compton scattering 

The mean free path of an electron due to the inverse Compton scattering, A/c, for 
the electron propagating through the isotropically distributed photon gas density 
n(w) is, 



P Jo 



(8) 



where the averaged over the angles cross section is 



(9) 



The total IC cross-section in the lab frame (see e.g. Ref. IT6| can be presented in 
the form 



07c 



3 or 

4x 



4 

X X" 



1 



1 



1 -^) ln(1 + X)+ 2 + x-2(l + x)2 



(10) 



where \ — 2wj(l — (3/j,). The exact solution of Eq.© can be presented in an analyt- 
ical form b u t ^he resulting expression contains a dilogarithm. A quite good 
approximation, with an accuracy better than 10% has been suggested in Ref. 1211 
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Fig. 3. The spectral energy distribution of upscattered photons (in arbitrary units) for 
different values of the parameter 6: 6 — » 0, b — 5, 10, 30, and b — > oo . 



T he sp ectral distribution of the scattered photons via IC for the isotropic photon 
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; x = — ; o = 4w7. 

7 7 



(12) 



The spectra of upscattered photons depend only on the parameter 6, as illustrated 
in Fig. [31 If the value of 6 is small, as in the Thomson regime, the upscattered 
photon population is predominantly in the low energy regime. If the value of 6 is 
large, corresponding to the Klein-Nishina regime, the population is shifted towards 
high energies. 



3. The method of calculations 

In order to study of formation of pair halos and to calculate the characteristics of 
their radiation, in this work we used the method of Monte Carlo simulations of 
interactions high energy electrons and photons with intergalactic radiation fields. 
This method provides an effective tool for study of development of electromagnetic 
cascades, and, in particular, allows precise calculations of both energy and spatial 
distributions of cascade particles for an arbitrary background photon field. The 
method of calculations presented below is described in Ref. [23j For the treatment 
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of interaction lengths and energy distributions of electrons and photons produced 
at the photon-photon and electron-photon interactions we used the corresponding 

JO Al 

total and differential cross-sections of these processes. The hit-or-miss method^and 
the rejection method^ have been used for simulations of the interaction lengths 
and the energy distributions of cascade particles. We assumed that the magnetic 
field in the intergalactic medium is sufficiently small, so the synchrotron cooling 
of electrons can be neglected. On the other hand, we assumed that the magnetic 
field is sufficiently large, thus the secondary electrons of all relevant energies are 
isotropised before they interact with ambient photons. These two conditions arc 
safely satisfied for a very broad range of the intergalactic magnetic fields between 
10- 7 G to 10- 12 G. 

The background photon fields used in our calculation consist of Cosmic Mi- 
crowave Background (CMB) and the Extragalactic Background Light (EBL). Both 
the CMB and EBL are functions of redshift. Since the background photon field at 
large cosmological redshifts significantly differ from that at the present epoch, the 
assumption of a consistent background during the development of the cascade can- 
not be held for the large redshift sources. On the other hand, the cascade simulation 
with evolving background photon field are computationally expensive. To make cal- 
culations faster, the simulation scheme was adjusted by simulating the cascade only 
at a certain epoch or in the limited region around a source of redshift z s , when 
one can assume constant background radiation with the EBL flux corresponding 
the cosmological epoch of z = z s . Obviously, the simulation region should be big 
enough to guarantee an adequate accuracy of calculations of the spatial distribu- 
tions of particles around the central source. At the same time the region should 
be relatively compact so the constant background approximation can be applied. 
For example, for a source at a distance of 1 Gpc, the gamma-ray production region 
around the source with an angular size 5° corresponds to the linear size rs 100 Mpc, 
so we can safely assume that the EBL is not noticeably changed across this size or 
over the corresponding time of < 3 x 10 8 yr. 

Beyond the "cascade region", we take into account the interactions with EBL 
only in terms of absorption (i.e. without of the cascade treatment). Thus the results 
presented in this paper are relevant only for the core of halos within angular size of a 
few degree. Although pair halos do extend beyond 5°, their cores, in fact, present the 
most interest, at least from the point view of detectability of these huge structures. 
Therefore in most of cases, when we are interested in the compact cores of halos, 
the calculations are done using this approximate procedure. Within this approach, 
to calculate the observed energy distribution dN/dE, we modify the distribution 
dN Zs jdE calculated in the " cascade region" at the epoch z s , by multiplying to the 
attenuation factor, exp[— rpp(e)]: 



(13) 
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Fig. 4. The attenuation factor exp(— rpp) for gamma photons emitted from sources at 
different redshifts as a function of energy for a specific EBL proposed by Primack et al. 
(2000). 



where Tpp is the PP optical depth, 



Here e z and u) z are the 7-ray and background photon energies at the redshift z, Hq 
is the Hubble constant; in this paper we use Hq= 65 km/s/Mpc. The attenuation 
factors at different redshifts are shown in Fig. 2J 

4. Results 

Below we present results of calculations of the spectral energy distributions (SEDs) 
and angular distributions of gamma-rays from pair halos around sources located at 
different redshifts. The main uncertainty of predicted energy spectra and morphol- 
ogy of pair halos is related to the uncertainties of EBL, especially in the mid and far 
infrared wavelengths. Here we use three different EBL models proposed by Malkan 
& Stecker (1998, hereafter M98j3, Primack et al. (2000, hereafter P00)Pand Kneiske 
et al. (2002, hereafter K02]P. One should note the recent developments related to 
the EBL at near infrared and optical wavelengths (see e.g. Ref. 8,9). However this 
wavelength band is not critical for formation of pair halos. On the other hand the 
above models represent quite broadly the possible range of EBL at most principal 
(from mid- to far- infrared) bands responsible for formation of pair halos around 
sources with redshifts up to z ~ 2. 

In Fig. [5] we present the Spectral Energy Distribution (SED) of a pair halo 
within different angles formed around a source with redshift z s = 0.129. This cor- 




(14) 



July 17, 2009 18:42 WSPC/INSTRUCTION FILE anant 



8 A. Eungwanichayapant and F.A. Aharonian 




2 3 4 2 3 4 

10 10 10 10 10 10 10 10 

ENERGY LGeV J ENERGY [GeV] 



Fig. 5. The SED of a pair halo at z = 0.129 within 6 < 10°, 1°, 0.3°, and 0.1° (left 
panel) and 0.1° < 9 < 0.3°, 0.3° < 9 < 1°, and 9 > 1° (right panel). In both panels 
the SEDs integrated over all angles are also shown. 

responds to the redshift of the blazar 1ES 1426+482, but is also representative for 
many TeV blazars reported by the HESS, MAGIC and VERITAS collaborations. 
For the EBL we adopted the POO model. The calculations are performed for a mo- 
noenergetic spectrum of primary gamma-rays with Eq = 100 TeV and luminosity 
L = 10 45 erg/s. 

The differential angular distribution of gamma-ray flux of the pair halo in dif- 
ferent energy bands is shown in Fig. [6] (left panel) . The angular distribution of the 
cumulative flux, J 7 (> 2 )/J 7 (all angles) is shown in Fig. [6] (right panel). The initial 
parameters are the same as in Fig. [5] As expected at higher energies = the halo 
looks more compact than at lower energies. This is clearly seen also in Fig. [Jj where 
the two-dimensional projections of pair halos at different energies are shown. 

4.1. The impact of the initial gamma-ray spectrum 

In order to study the impact of the initial gamma-ray spectrum on the basic char- 
acteristics of pair halos, we calculated the angular and spectral distributions of 
gamma-rays for a halo initiated by monoenergetic gamma-rays. The results are 
shown in Fig. [8] for three different energies of primary gamma-rays: Eo = 20, 100, 
and 500 TeV. One can see that at the intermediate (around the maximum of SED) 
and low energies both the energy and angular distributions are quite insensitive 
to the primary energy. The latter becomes important at higher energies, namely 
between 1 TeV and 10 TeV. In particular, with an increase of Eq the angular dis- 
tribution becomes sharper, especially below 0.1°. Also, with an increase of Eo, the 
energy spectrum beyond the maximum around a 100-300 GeV becomes significantly 
flatter. This tendency continues, however, until Eq ~ 100 TeV. At higher energies 
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Fig. 6. Left panel: differential angular distributions of the pair halo at z = 0.129 (in 
arbitrary units); Right panel: cumulative angular distributions at different energies. 



of primary gamma-rays, the results of the spectral and angular distributions of the 
halo radiation becomes almost insensitive to Eq. 



4.1.1. Power law sources 

In Fig. [9] we show the energy and angular distributions of of radiation of the pair 
halo for more realistic, power-law distribution of primary gamma-rays, Q(E) oc E~ T 
between 100 GeV and 100 TeV, and Q(E) = outside this interval. The results 
calculated for three different values of r=l, 1.5, and 2, but for the fixed total 
gamma-ray luminosity Lq = 10 45 erg/s are shown in Fig. 1101 The result are quite 
similar which is explained by the fact that for the chosen hard energy spectra of 
primary gamma-ray, the main contributions come from the energy interval close to 
£ max = 100 TeV. 

In Fig. [TO] we show the case of primary gamma-rays with fixed power-law photon 




E >10GeV E >100GeV E >500GeV E> 1 TeV 
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Fig. 7. Two dimensional map of a pair halo at z = 0.129 at different gamma-ray energies. 
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Fig. 9. Pair halos initiated by primary gamma-rays with power-law spectrum between 
-Emm = 100 GeV and E max = 100 TeV for different photon indices r=l, 1.5, and 2. Left 
Panel - SED; Right Panel - angular distribution. 



index T = 1.5, and for different maximum energies, E max = 10, 30, 100 TeV. It is 
interesting that the pair halo SEDs are more sensitive to -Emax between 10 to 30 
TeV, whereas the angular distributions are more sensitive to E max between 30 to 
100 TeV. 

4.2. Redshift Dependence 

When the redshift of the source exceeds 0.1 one has to take into account the evolu- 
tion of EBL in the past. The energy density of the CMB increases with z propor- 
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tional to (1 + z) 4 , while the peak of distribution shifts to higher energies by a factor 
of {l + z) (see FigfTTj). The EBL evolution is different from the CMB evolution since 
both the stellar-light and the dust component of EBL are produced at much later 
epochs with z < 3. 

In Fig. [11] we show the evolution of EBL, i.e. the EBL flux at different epochs z 
in accordance with the POO model (left panel) . The corresponding mean free paths 
are shown in the right panel of Fig. [11] One can see that at the epoch z = 2, the 
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Fig. 11. Left panel: EBL at different redshifts; Right panel: mean free paths of gamma-rays due 
to pair production (PP) and electrons due to the inverse Compton (IC) scattering at different 
epochs. 
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Fig. 12. Spectral luminosities of pair halos calculated for sources located at different red- 
shifts, z = 0.034, 0.129, 1 and 2, for the opening angle 1° (left panel). Angular distributions 
of gamma-rays at energies above 10 GeV (right panel). 

mean free paths of gamma-rays and electrons in the intergalactic medium were an 
order of magnitude shorter than at the present epoch. This should obviously result 
in significant shift of the SED towards low energies and more compact halos (in 
addition to the effects introduced by the redshift itself). This is demonstrated in 
Fig. 1121 where the spectral luminosities and angular distributions of radiation of 
pair halos are shown for different redshifts of central sources. The total luminosity 
of primary gamma-rays is assumed 10 erg/s. 

4.3. Dependence on the EBL models 

For energy spectrum of primary gamma-rays extended to E < 10 3 TeV, the e pairs 
are produced mainly through interactions with the EBL photons. Therefore, the 
choice of the EBL model is principal for calculations of characteristics of pair halos. 
To study the impact of the EBL on the pair halos, we calculated the parameters 
of pair halos for two different EBL models. The EBL from POO is based on semi- 
analytic calculations. All physical processes, including the formation and evolution 
of galaxies are included in the code. The EBL model M98 is based on empirical 
assumptions. It adopts an EBL flux at the present epoch (z = 0, and assumes time 
evolution proportional to (1 + z) 3A . These two EBL models for z = 0.129, as well 
as the corresponding mean free paths of gamma-rays are shown in Fig. [13] . 

In Fig. [14] we show the SEDs and angular distributions of a pair halo located at 
z = 0.129 calculated for two EBL models shown in Fig. [131 One can see that while 
the SEDs are quite similar at low energies, the energy spectrum corresponding to the 
POO model is significantly shifted compared to the SED predicted by the M98 EBL 
model. Also, the POO EBL model predicts sharper angular distribution, especially at 
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Fig. 13. Left panel: Spectral energy distributions of the POO and M98 EBL models at 
the present epoch, z — 0; Right panel: the mean free paths of gamma-rays for these two 
models, respectively. 



> 0.1°. This is explained by the higher EBL flux at mid and far- infrared wavelengths 
assumed in the POO model. 

For sources with large redshifts, z 3> 0.1, the energy and angular distributions 
of pair halos depends on the time-evolution of EBL. In fact, this dependence is 
stronger than the dependence on the EBL flux at the present epoch. To demonstrate 
this effect, in Fig. Qj}] we show the spectral and angular distributions of a halo at 
z = 2 calculated for two EBL models - POO and K02 '■ , which are characaterised 
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Fig. 15. The EBL fluxes at different redshifts predicted by the POO and K02 models. 

by significantly different cosmological evolution of the EBL flux. These two models 
predict relatively similar fluxes at present epoch (or at low redshifts, z < 0.1). 
However, at large cosmological epochs, z > 1, the difference of two models becomes 
as large as an order of magnitude. This is seen in Fig. [T5l where the EBL fluxes 
predicted by the POO and K02 model are shown for six different cosmological epochs. 
The corresponding SEDs and differential angular distributions are shown Fig. [16] 
and [lTj One can see that while at low redshifts two EBL models give relatively 
similar spectral and angular distributions of gamma-rays, at z — 2 the difference in 
predictions of two EBL models becomes significant. 
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Fig. 16. Spectral luminosity distributions of the pair halos corresponding to each EBL 
models 

5. Summary 

Interactions of high energy gamma-rays with EBL initiate electroimagnctic cascades 
which lead to formation of electron-positron halos around powerful nonthermal ex- 
tragalactic sources. The pair halos contain unique information about the flux of 
EBL and its time-evolution in the past. We studied the dependence of radiation of 
pair halos on several key parameters, in particular on the energy spectrum of the 
primary gamma-rays, the source redshift, and the EBL flux. The results presented 
in this paper can be used for search for these giant structures which can serve as 



July 17, 2009 18:42 WSPC/INSTRUCTION FILE anant 



16 A. Eungwanichayapant and F.A. Aharonian 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

Q 2 [degree 2 ] 



Fig. 17. 



unique cosmological candles. 
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